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LONG-TERM  GOALS 

To  improve  our  undestanding  and  predictive  capability  for  three  dimensional,  time -dependent  ocean 
circulation. 
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OBJECTIVES 


Our  main  objective  is  to  develop  modeling  approaches,  theory  and  observational  diagnostics  that  will 
allow  a  better  understanding  of  the  role  of  processes  that  are  intrinsically  three  dimensional  and  time 
varying.  We  hope  to  use  this  knowledge  to  inform  parameterizations  and  predictions,  and  also  to  guide 
observational  strategies,  including  strategies  for  autonomous  underwater  sensors  and  vehicles. 

APPROACH 

Our  approach  centers  on  the  Lagrangian  view  of  ocean  circulation;  that  is,  the  study  of  fluid  motion 
through  tracking  of  fluid  trajectories.  The  methods  of  analysis  are  based  in  the  mathematical  field  of 
dynamical  systems. 

WORK  COMPLETED 

We  have  developed  a  number  of  methods  for  identifying  and  computing  internal  boundaries  and  other 
coherent  features  that  drift  and  evolve  with  ocean  currents  and  that  separate  qualitatively  different 
regions.  These  structures  provide  a  template  for  understanding  fluid  pathways,  mixing,  and  transport 
processes  in  3D.  These  methods  have  been  tested  using  a  variety  of  models  ranging  from  idealized  to 
realistic.  We  have  also  developed  diagnostics  that  indicate  when  two-dimensional  analysis  suffices. 
Another  major  thrust  involves  the  development  of  methods  for  assimilating  Lagranigan  data  into  ocean 
models. 

RESULTS 

We  note  that  our  program  was  the  subject  of  a  3-year  review  at  ONR  in  May  of  this  year.  All 
presentations  from  that  meeting  can  be  downloaded  from  the  our  MURI  program  website: 

http://www.whoi.edu/ocean3dplusl/ 

under  MURI  meetings. 

New  methods  for  computing  barriers  and  other  Lagrangian  coherent  structures  (LCS)  include  the 
ergodicity  defect  (a  measure  of  trajectory  complexity),  Koopman  Operator  techniques  (which  rely  on 
properties  measured  along  trajectories),  techniques  that  identify  regions  of  high  hyperbolicity 
(corresponding  to  rapid  fluid  stretching  and  folding),  and  methods  that  calculate  the  rate  of  separation 
of  fluid  from  it  neighboors.  Examples  are  depicted  in  Figures  1-5.  All  have  been  proved  capable  of 
identifing  key  Lagrangian  structures  in  the  flow  field  and  the  choice  of  a  particular  method  will  depend 
on  the  nature  of  the  model  or  observation  information  available  We  have  also  made  progress  in 
identifying  the  kind  of  barriers  that  can  exist  in  fully  3D,  time -dependent  flows  (Figures  5-7  and  9). 
Our  serarch  for  barriers  is  being  extended  to  more  realistic  3D+1  simulations  (various  figures). 
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Figure  1.  A  method  for  calculation  of  LCS  using  measures  of  trajectory  complexity: 
work  by  Rypina,  Pratt,  Scott,  and  Brown. 


Gulf  Oil  Spill  Prediction 


Hypwgraph  0  Observed  oil  spread 


(A)  Ocean  hypergraph  map  in  front  of  the  Bioa-Pensacola  shoreline  on  25  Way.  forecasting 
strong  oil  incursion  toward  the  coastline  (circled)  in  the  fd.  lowing  3  days,  iflj  NONA'S  oil  spread  estimate 
in  front  of  the  Biloxi-Pensacola  shoreline  on  27  Ma  y.  Tt»  major  directions  of  oil  spread  were  predicted  hy 
;he  hypergraph  map  2  day*  earlier.  Ih#  ol  reached  die  shore  several  days  later,  on  2  June.  (C)  Ocean 
hypejgraph  map  in  I  rent  of  Pensacola  on  G  June.,  forecasting  a  strong  oil  mixing  event  in  fro  nt  of  Ihe 
shoreline  and  extension  of  the  oil  slick  toward  Panama  City  Beach  in  the  following  3  days.  ID)  NOAA’s  oil 
spread  estimate  on  10  June  in  front  of  Pensacola.  The  oil  developed  a  large  slick  forecasted  by  the 
hypengraph  map  2  days  eadier  and  continued  to  flow  toward  Panama  City  Beach. 
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Figure  2.  Mesohyperbolicity  identifies  areas  of  strong  filamentation  (Mezic). 
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Hyperbolic  structure  in 
temperature  and  flow 
fields  during  Loop  Current 
Ring  pinchoff 

Mixing  boundaries  (FTLEs)  over 
0-200  meters  computed  from 
GOM  HYCOM  model  velocities 
shown  as  red  and  blue  curves. 

Mixing  boundary  geometry  is 
nearly  vertical  over  this  depth 
range. 

Saddle  points  in  the  temperature 
field  are  found  near  the 
stagnation  point  (and  FTLE 
intersection  point)  at  all  depths. 


June  10,  2012 


O  Stagnation  points 


|  ,  ,  ,  _3  ftwHttiMitifcfHa  O  Temperature  saddle  points 

•OS  03  DjJ  0.7 


Sulman  et  al.  Hyperbolicity  in  temperature  and  flow  fields  during  the  formation 
of  a  Loop  Current  ring,  Nonlinear  Proc.  Geophys,  in  press,  2013. 


Figure  3.  Use  of  mesohyperbolicity  to  identify  barriers  and  strongly  stirred  regions. 


MeziC  Research  Group 


Computational  Geometric  Analysis 

Goal:  detect  locally  connected  subsets  of  EQ  to  identify  trajectories  that  lie 
within  a  regions  with  local  integrals  of  motion. 


Trajectories 


Ergodic  Quotient 
(extrinsic  coordinates) 


Ergodic  Quotient 
(intrinsic  coordinates) 


Coherent 

Structures 


Compute  numerical  diffusion  modes  on  the  ergodic  quotient  to 
obtain  low-dimensional  representation  of  the  ergodic  quotient. 

Geometric  clustering,  e.g.,  k- means,  in  diffusion  space  can  then  extract  connected 
subsets.  Color  trajectories  based  on  the  membership  in  clusters. 


Figure  4.  The  Ergodic  Quotient  is  based  on  Lagrangian  averages  of  scalar  quantities. 
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Figure  5.  Comparison  with  known  examples  such  as  ABC  flow  is  favorable. 


Figure  6.  Idealized  models  of  fully  3D  eddy  circulations  tell  us  what  Lagrangian  barriers  look  like  in 
3D.  In  this  case,  the  barriers  consist  of  tori,  including  twisted  hula  hoops.  Sandwiched  between 
the  surfaces  are  bands  and  islands  of  chaos.  (From  Pratt,  Rypina,  Ozgokmen,  Wang, 

and  Childs.  JFM  2013. 
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Object  resulting  from  double  resonance 


Figure  7:  When  the  flow  is  fully  3D  and  time-dependent,  the  barriers  become  exotic. 
In  this  case  the  cocoon-like  structure  pulsates  and  rotates  in  the  eddy. 

Do  the  3D  FTLEs  detect  any  coherent  features  in  a  highly  turbulent  situation? 


Figure  8:  The  methods  used  to  compute  LCS  also  provides  information  about  changes  in  the 
stirring  distribution  in  the  cross  section  of  a  fully  3D  turbulent  rising  oil  plume.  (Ozgokmen) 
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Our  team  is  also  exploring  the  application  of  LCS  methodology  in  connection  with  other  realistic 
models.  This  work  emphsizes  the  role  of  the  submesoscale  (Figs.  10-14),  which  is  commonly 
unresolved  in  ocean  models. 


Figure  9.  We  can  get  a  path  of  understanding  into  fully  3D,  time  dependent  systems  but  looking  a 
weakly  3D  systems  for  long  time  intervals.  (Llewellyn  Smith) 

3D  FTLE: 


Submesoscale  phase: 


Mesoscale  phase: 


Clearly  different  turbulent  coherent  structures: 

shallow  submesoscale  eddies  vs  deep  mesoscale  features... 


Figure  10.  We  are  exploring  the  effects  of  submesoscale  motions  (left),  which  are  unresolved  in 

most  models.  (Ozgokmen) 
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Turbulent  exchange  quantified  from  tracer  release: 


time  (days) 


Figure  11.  An  order  of  magnitude  increase  in  frontal  turbulent  diffusivity  after  the  deep  instability 
kicks  in;  submesoscale  overcome  by  mesoscale  in  this  case. 


Cohabitation  of  Mesoscale  and  Submesoscale: 


NQAA  m02  AVHRR  SST  (C).  01  *F#b  2QI  2  0 1  56  00  UT( 


■76  -74  -72  *70  -66 

LonE 


Figure  12.  The  instabilities  that  cause  star  eddies  reside  in  the  submesoscale  (Ozgokmen). 
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3D  FTLEs  From  LES  of  Oceanic  Eddies: 


Figure  13.  We  are  searching  for  barriers  in  realistic  3D,  time-dependent  eddies.  Do  we  see 
generalizations  of  the  tori  found  in  idealized  models?  (Ozgokmen,  Pratt,  Rypina,  and  Wang.) 


Which  Scales  Control  Transport  in  the  Ocean? 


Poje  et  al.  (2010) 
Haza  et  al.  (2012) 


/I  2  degree  HYCOM 


Image  credit:  Haza 


/48  degree  HYCOM 


(I)  Are  the  long-living,  slow  mesoscale  features  enough  to  compute  transport? 

(II)  Or,  rapidly-evolving,  smaller  submesoscale  transport  barriers  are  needed? 

Figure  14.  Comparison  between  LCS  are  different  model  resolution.  The  lower  panel 

resolves  the  submesoscale. 
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In  the  area  of  Lagrangian  Data  Assimilation  (Fig.  15),  we  are  developing  methods  that  are  more 
efficient  and  accurate,  and  that  work  in  3D.  This  work  attempts  to  use  information  obtained  from  LCS 
analysis  in  order  to  determine  where  to  best  lauch  floats  and  drifters  and  to  optmize  navigation  strategy 
for  autonomous  underwater  vehicles  Fig.  16. 


Augmented  system 


X  = 


Append  equations  for  drifters  (floats,  gliders,  AUVs) 

Ide,  Jones  and  Kuznetsov  (2002) 

—  augmented  state  vector 


X  ' 


D 


dx 


dt 

dx\ 


—  =  M  (x  ,J)  —  flow  equations 


D 


dt 


=^(*0.4.0 


advcction  equation 


Apply  filtering  to  augmented  system: 

1.  Ensemble  Kalman  Filter 

2.  MCMC  (Particle  Filter) 


Figure  15.  The  general  ideal  behind  Lagrangian  data  assimilation  work  of  Jones, 

Spiller  and  colleagues. 
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Instrument  control  to  optimize  information 


Add  simple  control: 

/(z)  =  (f,0) 

0,5 
0,4 
0,3 
0.2 
0,1 

0.0 

0 

yi  -  x(tk)  4-  m  ~  Ar(0,  a-2/);  Jfe  =  I,  * . K/2  no  control 
y*  =  x(4)  +  fjfc,  Tjjt  ~  N{0,  a3/),  t  -  K/2,  ,,.,K  control 

Want:  P(v\y\y2) 


dv 

dt 

dz 

dt 


=  0 


—  =  v(z)  +  f(z) 


Figure  16.  Is  it  possible  to  optimize  navigation  strategies  for  AUV’s  in  order  to  collect  the 
best  data  for  LDA  ?  (Jones,  Spiller  and  colleagues.) 


IMPACT/APPLICATIONS 

The  Lagrangian  view  offers  insights  into  transport  and  mixing  processes  that  are  different  from  those 
gained  in  the  Eulerian  frame.  These  insights  naturally  lead  to  better  physics-based  parameterizations, 
such  as  turbulent  eddy  coefficents.  (One  of  the  criticisms  of  Eulerian-based  parameterizations  is  that 
they  depend  only  on  the  local  properties  of  the  flow  and  do  not  consider  patches  of  turbulence  that  are 
advected  into  the  local.)  Lagrangian  methods  can  also  identify  time-dependent  barriers  that  may  be 
important  in  predicting  regions  that  polutants,  including  oil  and  radioactive  material,  have  difficulty 
crossing. 

The  Lagrangian  approach  to  data  assimilation  will,  in  principle,  allow  information  that  has  been 
previoiusly  discarded  to  be  used  in  prediction  and  state  estimation.  In  particular,  previous  assimilation 
schemes  that  use  sequences  velocity  measrments  from  drifting  sensors  do  not  assimilate  the 
information  that  each  measurement  comes  from  the  same  trajectory. 

TRANSITIONS 

The  tools  being  developed  will  aid  in  the  prediction  of  catestrophic  release  of  oil,  radioactive  material 
or  floating  debree  into  the  ocean. 


12 


RELATED  PROJECTS 


A  British  Petrolium-funded  project  on  the  Gulf  of  Mexico  oil  spill.  Prof.  Tamay  Ozgokmen  is 

involved  along  with  several  other  Pis. 

PUBLICATIONS 

Amey,  C.,  B.  Sundaram,  and  A.  C.  Poje.  Mixing  in  mixed  phase  spaces.  To  be  published  in  Physics  of 
Fluids,  May  2013. 

Apte,  A.,  and  C.K.R.T.  Jones,  2013.  The  impact  of  nonlinearity  on  Lagrangian  data  assimilation. 
Nonlinear  Processes  in  Geophysics,  20,  329-341. 

Apte,  A.,  E.  Spiller,  and  C.K.R.T.  Jones.  Assimilating  en-route  Lagrangian  observations.  Tellus  A 
[accepted,  refereed], 

Balasuriya,  S.,  G.  Froyland,  and  N.  Santitissadeekom,  2014.  Absolute  flux  optimising  curves  of  flows 
on  a  surface.  Journal  of  Mathematical  Analysis  and  Applications,  409,  1 19-139. 

Budisic,  Marko,  and  Igor  Mezic,  2012.  Geometry  of  the  ergodic  quotient  reveals  coherent  structures  in 
flows.  Physica  D,  241(15),  1255-1269,  10.1016/j.physd.2012.04.006. 

Chabreyrie,  R.,  and  S.  G.  Llewellyn  Smith.  Short-  and  Long-Time  Structures  in  a  Three-Dimensional 
Time  Dependent  Flow.  Submitted  to  Physics  of  Fluids. 

Cox,  G.  Large -time  uniqueness  in  a  data  assimilation  problem  for  Burgers'  equation,  2012.  Inverse 
Problems,  http://arxiv.org/abs/1207.4782vl. 

Cox,  G.  Uniqueness  and  non-uniqueness  of  minimizers  in  variational  data  assimilation.  Submitted  to 
Physica  D. 

Cox,  G.  The  L  essential  spectrum  of  the  2D  Euler  operator.  Submitted  to  Journal  of  Mathematical 
Fluid  Mechanics. 

Froyland,  G.,  C.  Horenkamp,  V.  Rossi,  N.  Santitissadeekom,  and  A.  Sen  Gupta,  2012.  Three- 

dimensional  characterization  and  tracking  of  an  Agulhas  Ring.  Ocean  Modelling,  52-53,  pg.  69- 
75. 

Griffa,  A.,  A.  C.  Haza,  T.  M.  Ozgokmen,  A.  Molcard,  V.  Taillandier,  K.  Schroeder,  Y.  Chang,  and 
P.M.  Poulain,  2013.  Investigating  transport  pathways  in  the  ocean.  Deep  Sea  Research  II,  85,  81- 
95. 

Haza,  A.  C.,  T.  M.  Ozgokmen,  A.  Griffa,  Z.  D.  Garraffo,  and  L.  Piterbarg,  2012.  Parameterization  of 
particle  transport  at  submesoscales  in  the  Gulf  Stream  region  using  Lagrangian  subgridscale 
models.  Ocean  Modelling,  42,  31-49. 

Huntley,  H.  S.,  B.  L.  Lipphardt,  Jr.,  and  A.  D.  Kirwan,  Jr.,  2011.  Surface  Drift  Predictions  of  the 
Deepwater  Horizon  Spill:  The  Lagrangian  Perspective.  Geophysical  Monograph  195: 
Monitoring  and  Modeling  the  Deepwater  Oil  Spill:  A  Record-Breaking  Enterprise,  pg.  179-195. 

Loire,  S.,  P.  Kauffmann,  I.  Mezic  and  C.  D.  Meinhart,  2012.  A  theoretical  and  experimental  study  of 
AC  electrothermal  flows.  J.  Phys.  D:  Appl.  Phys.,  45,  185301. 

Mensa,  J.,  A.  Griffa,  Z.  Garraffo,  T.M.  Ozgokmen  and  M.  Veneziani:  Seasonality  of  the  submesoscale 
dynamics  in  the  Gulf  Stream  region.  Ocean  Dynamics,  submitted:  March  2013. 


13 


Ozgdkmen,  T.M.,  A.C.  Poje,  P.F.  Fischer,  H.  Childs,  H.  Krishnan,  C.  Garth,  A.  Haza,  and  E.  Ryan, 

2012.  On  multi-scale  dispersion  under  the  influence  of  surface  mixed  layer  instabilities  and  deep 
flows.  Ocean  Modelling,  56,  16-30  DOI:  10.1016/j.ocemod.2012.07.004. 

Ozgokmen,  T.  M.,  A.  C.  Poje,  P.  F.  Fischer,  and  A.  C.  Haza.  Large  eddy  simulations  of  mixed  layer 
instabilities  and  sampling  strategies,  2011.  Ocean  Modelling,  39(3-4),  311-331,  DOI: 

1 0. 1 0 1 6/j  .ocemod.20 11.05.006. 

Piterbarg,  L.I.,  V.  Taillandier,  and  A.  Griffa.  Investigating  frontal  variability  from  repeated  glider 
transects  in  the  Ligurian  Current  (North  West  Mediterranean  Sea).  J.  Mar.  Sys.,  submitted  (April 
2,2013). 

Pratt,  L.  J.,  I.  Rypina,  T.  M.  Ozgokmen,  and  Y.  Bebieva.  Chaotic  advection  in  a  rotating-cylinder 
flow.  Submitted  to  J.  Fluid  Mech. 

Pratt,  L.,  I.  Rypina,  T.M.  Ozgokmen,  H.  Childs,  and  Y.  Bebieva.  Chaotic  advection  in  a  steady,  three- 
dimensional,  Ekman-driven  eddy.  Journal  of  Fluid  Mechanics,  in  revision,  March  2013. 

Rypina,  I.,  S.  E.  Scott,  L.  J.  Pratt,  and  M.  G.  Brown,  2011.  Investigating  the  connection  between 
complexity  of  isolated  trajectories  and  Lagrangian  coherent  structures.  Nonlin.  Processes 
Geophys.,  18,  977-987. 

Rypina,  I.I,  L.  Pratt,  and  I.  Mezic.  Resonance  phenomena  in  3D  time-dependent  volume-preserving 
flows  with  symmetries.  In  preparation  to  Nonlin.  Processes  Geophys.,  Special  Issue  on 
"Nonlinear  Processes  in  Oceanic  and  Atmospheric  Flows". 

Santitissadeekom,  N.,  E.T.  Spiller,  C.K.R.T.  Jones,  R.  Rutarindwa,  L.  Liu,  and  K.  Ide.  Observing 
System  Simulation  Experiments  of  Cross-Layer  Lagrangian  Data  Assimilation.  Submitted  to 
Dynamics  of  Atmospheres  and  Oceans,  January,  2013. 

Schroeder,  K.,  J.  Chiggiato,  A.  C.  Haza,  A.  Griffa,  T.  M.  Ozgokmen,  P.  Zanasca,  A.  Molcard,  M. 
Borghini,  P.M.  Poulain,  R.  Gerin,  Z.  Zambianchi,  P.  Falco,  and  C.  Trees,  2012.  Targeted 
Lagrangian  sampling  of  submesoscale  dispersion  at  a  coastal  frontal  zone.  Geophys.  Res.  Lett., 

39,  LI  1608,  doi.10.1029/2012GL051879. 

Sulman,  M.  H.  M.,  H.  S.  Huntley,  B.  L.  Lipphardt,  Jr.,  and  A.  D.  Kirwan,  Jr,  2012.  Out  of  Flatland: 
Three-Dimensional  Aspects  of  Lagrangian  Transport  in  Geophysical  Fluids.  Geophysical 
Monograph  200:  Lagrangian  Modeling  of  the  Atmosphere,  p.  77-84. 

Sulman,  M.  H.  M.,  H.  S.  Huntley,  B.  L.  Lipphardt  Jr.,  and  A.  D.  Kirwan  Jr.,  2013.  Leaving  Flatland: 
Diagnostics  for  Lagrangian  coherent  structures  in  three-dimensional  flows.  Physica  D,  258:  77- 
92. 

Sulman,  M.  H.  M.,  H.  S.  Huntley,  B.  L.  Lipphardt  Jr.,  G.  Jacobs,  P.  Hogan,  and  A.  D.  Kirwan  Jr., 

2013.  Hyperbolicity  in  temperature  and  flow  fields  during  the  formation  of  a  Loop  Current  ring. 
Nonlinear  Proc.  Geophys.  [in  press] 

Valentine,  D.  A.,  I.  Mezic,  Senka  Macesic,  Nelida  Cmjaric-Zic,  Stefan  Ivic,  P.  J.  Hogan,  V.  A.  doi: 

10.1 073/pnas.  1 1 08820 1 09  water  hydrocarbon  irruption.  Proceedings  of  the  National  Academy  of 
Sciences,  doi:  10. 1073/pnas.l  108820109. 


14 


